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Abstract 

The space-time conservation element and so- 
lution element (CE/SE) method is employed to 
numerically study the near-field axisymmetric 
screech-tone noise of a typical underexpanded 
circular jet issuing from a sonic nozzle. For 
the computed cases, corresponding to fully ex- 
panded Mach numbers of 1 . 10, 1.15 and 1.19, 
the self-sustained feedback loop is automat- 
ically established. The computed shock-cell 
structure, acoustic wave length, screech tone 
frequencies, and sound pressure levels are in 
good agreement with experimental results. 


1 Introduction 

An underexpanded supersonic jet radiates mixing 
noise, broadband shock-associated noise, as well as 
screech tones under certain conditions. These compli- 
cated and technologically important physical phenomena 
have been the topic of many experimental and theoretical 
investigations, see the review papers [ 1-4] for a compre- 
hensive list of references and discussion. Simply put, the 
mixing noise is directly associated with large-scale struc- 
tures, or instability waves, in the jet shear layer; whereas, 
the broadband shock-associated noise and screech tones 
are associated with the interaction of these waves with 
the shock-cell structure in the jet core. The screech tones 
arise due to a feedback loop, Le. % part of the acous- 
tic waves generated by the wave/shock-cell interaction 
propagate upstream and re-generate the instability waves 
at, or in the vicinity of, the nozzle lip. Screech tones are 
of particular interest not only because of general noise- 
reduction concerns, but also because of potentially de- 
structive structural interaction(s) leading to sonic fatigue. 
The feedback loop leading to screech tones is sensitive 
to small changes in the system conditions, and the un- 
derstanding of the phenomena is to date mostly based on 
experimental observations [5-9]. 


Reliable direct numerical simulation of jet screech 
noise has up to quite recently not been feasible. The 
numerical scheme is required, on one hand, to resolve 
the acoustic waves without too much dispersion error 
and numerical dissipation, while on the other hand, it is 
required to capture shocks, or other nonlinear phenom- 
ena, near or inside the jet correctly. In addition, non- 
reflecting boundary conditions must be implemented, 
which is more difficult to accomplish in the near field 
than in the far field. 

Shen and Tam [10,1 1] obtained excellent results in a 
direct numerical simulation of screech for circular jets, 
using the well-known DRP scheme. A medium size com- 
putational domain of 35 D long in the streamwise direc- 
tion and HD in radial direction is employed in their 
computation, where D is the nozzle diameter. 

The space-time conservation element and solution el- 
ement method (CE/SE) [12-15] is another scheme that 
meets the above requirements. As demonstrated in 
our previous papers, the CE/SE scheme is well suited 
for computing waves on compressible shear flows [16] 
as well as vorticity/shock interactions [17], both being 
cornerstones of the jet-noise phenomena. Because of 
the CE/SE non-reflecting boundary conditions (NRBC), 
which are based on dux balance , a much smaller near 
field computational domain (12D x 8/7) can be used in 
the present numerical simulation. 

The unstructured-grid axisymmetric CE/SE scheme 
used here is briefly discussed in Section 2. Section 3 
illustrates the initial and boundary conditions for the 
screech problem, in particular, the CE/SE NRBC. The 
numerical results are presented and compared to avail- 
able experimental findings [5-9] in Section 4. Conclu- 
sions are drawn in Section 5. 

2 The Unstructured Axisymmetric CE/SE 
Navier-Stokes Solver 

In this Section, several new developments of the 
CE/SE numerical method are applied: (a) use of an 
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unstructured grid, which helps to facilitate the numeri- 
cal implementation [15]; (b) the Navier-Stokes CE/SE 
solver; and (c) the axisymmetric LES (large eddy sim- 
ulation), which accounts for the strong momentum ex- 
change between the jet and the ambient flow and results 
in a much improved numerical simulation. These new 
developments are described in the following subsections. 


0,3 = (7 - 1 )U 4 + [(3 - 7 )Ui - (7 - 1)E£] /2 U u 
0,4 = 'tUsU.IU, - (7 - 1)03 [02 + 02] /20?, 

and the viscous fluxes are: 

F v i = 0, F v2 = /x(2u a: - §V-V), 


2.1 Conservation Form of the Unsteady 
Axisymmetric Navier-Stokes Equations 

In general, the CE/SE method systematically solves a 
set of integral equations derived directly from the physi- 
cal conservation laws and naturally captures shocks and 
other discontinuities in the flow. Both dependent vari- 
ables and their derivatives are computed simultaneously. 
As a consequence, the flow vorticity can be obtained 
without reduction in accuracy. Non-reflecting bound- 
ary conditions (NRBCs) are also easily implemented be- 
cause of the flux-conservation formulation. 

Consider a dimensionless conservation form of the un- 
steady axisymmetric Navier-Stokes equations of a per- 
fect gas. Let p, u, v, p, and 7 be the density, streamwise 
velocity component, radial velocity component, static 
pressure, and constant specific heat ratio, respectively. 
The axisymmetric Navier-Stokes equations then can be 
written in the following vector form: 

Ut + F x + G y = Q, (1) 


Fv3 = MTr + Uy), 

F v 4 - /x[2 uu x 4- (Uy 4- v x )v - |(V-V)u+ 

_ u<2 + )} 

Prdy^Ux 2 jJ ’ 

G v y = 0, G v 2 = p(v x 4- u y ), 

G v3 = p(2v y - f V.V), 

G V 4 = p[2 Wy 4- (u y 4- v x )u “ §(V*V)i?4- 

7 d U4 u 2 4- v 2 . 1 

~P~rdy^U[ 2 

where u, r, u y , are respectively the x— and y— 
flow velocity components and their derivatives, which 
can be written in terms of the conservative variables 
Ui,U 2 ,Us and t/ 4 , with Pr (= 0.72) being the Prandtl 
number, p the viscosity, the velocity divergence 


where :r, y ^ 0, and t are the streamwise and radial co- 
ordinates and time, respectively. The conservative flow 
variable vector U and the flux vectors in the streamwise 
and radial directions, F and G, are given by: 


U 



F = 


f Fi \ 

F-2 

F \ 

\fJ 


G = 



with 

U\ = p, U2 = pu , Us — pv, 

04 = p/(7 — 1) + p(u 2 + v 2 )/2. 

The flux vectors are further split into inviscid and viscous 
fluxes: 

F — Fj — F V ,G — Gj — G v , 

where the inviscid fluxes are the same as in the Euler 
equations: 

Fn - 02 , 

0,2 = (7 - 1)04 + [(3 - 7 )02 - (7 - 1)0?] /20i, 
F i3 = 0203/01, 

0,4 = 70204/01 - (7 - 1)02 [ 0 ? + 0 ?] / 2 Ul 
Gil = 03 , 0,2 = 0203 / 01 , 


V-V = u x -f Vy 4- vfy. 

The right hand source term Q is the same as in the 
axisymmetric Euler equations [18]: 


Q 



where 


Q 1 = -Us/y, Q2 — — ^Us/Uiy, 


Qs = -0|/0iy, <? 4 = -Gi/y. 

By considering ( x , t/, t) as coordinates of a three- 
dimensional Euclidean space, £3, and using Gauss’ di- 
vergence theorem, it follows that Eq. ( 1 ) is equivalent to 
the following integral conservation law: 



■dS = 



rn — 1, 2, 3,4, (2) 


where 5(V r ) denotes the surface around a volume V in 

£3 and H m — (F m , G m , U m ) . 
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a hexagon cyfcider AD6ECFA - A’D’B’EC’FA 
w»>te3CEs:QADB-0’AD’B’, OBEC-O’B’E’C , 
andOCFA-aCFA’ 



b. projection of E3 in the x-y plane 

Figure 1 : CE/SE unstructured grid 


2.2 Unstructured Grid for CE/SE 

The CE/SE scheme is constructed to take advantage of 
an unstructured triangle grid. The unstructured geome- 
try used with the CE/SE scheme is illustrated in Fig. 1. 
Here, AABC is a typical triangle cell and D,E,F are 
the triangle centers of the neighboring cells. The flow 
variables at the previous time step are stored at the tri- 
angle cell centers. Three quadrilateral cylinders (conser- 
vation elements) are formed by the edges that connect 
the vertices and the center of the triangle and its three 
neighbors. The projection of these cylinders in the x-y 
plane are shown in Fig. 1(b) and the CE's are shown in 
Fig. 1(a). In the space-time £3 space, Eq. (2) is applied 
to the hexagon cylinder (of volume V) that consists of 
these 3 quadrilateral cylinder CE s, see Fig 1(a). The 
discrete approximation of Eq. (2) is then 



• dS - V(Q m ) 


n + 1 

j 


(3) 


form = 1,2, 3. 4, where H* n = (F* n , G* n , U* n ). The 
right-hand side of Eq. (3), in general, is the volume V 
times the source term evaluated at an appropriate Gaus- 
sian quadrature node. Here, the Gaussian quadrature 
node is the center of the hexagon A' D f B' E'C f F f at the 
time level n + 1, /. e . , O'. 

In the CE/SE scheme, the above flux conservation re- 
lation, Eq. (3), in space-time is the only mechanism that 
transfers information between node points. A conserva- 
tion element CE (here, quadrilateral cylinders) is the fi- 
nite volume to which Eq. (3) is applied. Discontinuities 
are allowed to occur in a conservation element. A so- 
lution element SE associated with a grid node ( e.g ., D, 
E , or F in Fig. 1 ) is here a set of interface planes in £3 
that passes through this node (e.g. DAA'D\ DBB f D\ 
EBB'E\ ECC'E\ etc.). Each surface S(CE) is made 
up of segments belonging to two neighboring CE' s. 
Within a given solution element SE(j , n), where j, n are 
the node index and time step, respectively, the flow vari- 
ables are not only considered continuous but are also ap- 
proximated by the linear Taylor series expansions: 


U*(x, y, t: j, n) = U r j -f (U x Yj(.r - Xj)+ 

(U y )Ky- yj ) + (U t )?(t-t n )< ( 4 ) 

F*(x.y,tj,n) = F] + (F r )?(x - xj)+ 

(F y )](v-yj) + (F t )y(t-t n ). (5) 

G*(.r, yj j. n) = G’j + (G x )j(.r -- .ty) + 

(Gj;(y~ yj )F(G t YJ(t-t n l (6) 

where j is the node index of D , E or F. The partial 
derivatives of F and G can be related to the correspond- 
ing one of U by using the chain rule with U t provided 
by Eq. (I ). The surface flux can be calculated accurately 
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and easily by first evaluating the flux vectors at the ge- 
ometrical center of the surface through the above Taylor 
series expansions. 

At time level rc, the solution variables U , U x , and 
U y are given at the three nodes D . E , F in Fig. 1 and 
U , U x and U y at O' at the new time level n + 1 are 
to be computed. In principle, each of the 3 CE s pro- 
vides 4 scalar equations when Eq. (3) is applied to the el- 
ement. Hence, the 12 scalar equations needed for the 12 
scalar unknowns at O' are available. All the unknowns 
are computed based on these relations. No extrapola- 
tions (interpolations) across a stencil of cells are needed. 
In practice, however, the dependent variables are first ob- 
tained at the hexagon center at the new time level n -f 1 
using Eq. (3) and then their values at the center O' of the 
triangle ABC can be obtained by a linear Taylor series 
expansion. Details can be found in [15]. 

With an unstructured grid, the CE/SE procedure is 
simplified and more easily adapted to complicated ge- 
ometries. Only a single mesh is needed as compared 
to the spatially staggered mesh used in the structured- 
grid case and the time-marching is completed in one step 
rather than two. Furthermore, the simple non-reflecting 
boundary conditions described previously [16-19] still 
work well with an unstructured grid. More details about 
the unstructured CE/SE method can be found in [15]. 
The weighted a - e CE/SE scheme is used here. 

2.3 Treatment of the Source Term 

The treatment is identical to the one used in [18] and is 
briefly reiterated here. Since the source term Q itself is a 
function of the unknown U, a local iterative procedure is 
needed to determine U . The discretized integral equation 
(3) reduces to the form 

U -Q(U)At = U H , (7) 

where U // is the local homogeneous solution ( Q = 0 
locally). Note that U h only depends on the solution at 
the previous time step, i.e ., U // is obtained using explicit 
formulas. A Newton iterative procedure to determine U 
is then 

u il+ 11 = u {l) - (|i)- l [#(C/ (<) ) - U„l 

where i is the iteration number and 


<!>(£/) — U~ Q(U)At. 


Normally, U at the previous time step is a good initial 
guess U i0) and the procedure takes about 2-3 iterations 
to converge. The Jacobian matrix is given by 


OU 


( 

\ 



0 

1 + LLi&l 
1 + ihv 

0 


A 2 


At 

y 

U 2 At 
thy 


l _l_ LhAt 

ihv 

A; 


0 


1 + ) 

2f r i y / 


where 


A\ 


Us At (7-1 m + Ui) 

— — 71/4 — 

U iy U 1 


A-2 = -(7 - 1 ) 


U 2 U 3 At 

U\V ’ 


4 A* r 7-1^1 +3t/2 7l7 4l 

A3 = T [ "^ U^ + Ih ] 


The inverse of the Jacobian, i.e., ( Jjj) -1 can easily be 
derived analytically for this particular case, thus, leading 
to a savings in CPU time. 


2.4 Large Eddy Simulation (LES) 

A simplified LES procedure similar to those used in [20] 
is adopted here to account for the strong momentum ex- 
change in the shear layer. In the LES, a simple Smagorin- 
sky subgrid scale model is used for the eddy viscosity: 


th = (C s A) 2 (2S l3 S tj ) l/2 , 


where 


— 9 ( 


duj du 1 


+ 


dx , 


). 


2 v dx j 

with A — (AxAtj)^ 2 and C s = 0.1. Then, // + fit 
replaces fi in the actual computation. In practice, the 
LES model was only applied for radii less than 1.5 D and 
4 D downstream of the nozzle exit plane. 


3 The Jet Screech Noise Problem 

Consider a circular jet as sketched in Fig. 2. The flow 
at the start of the computations is choked, i.e., nozzle exit 
Mach number, M e , is unity, and the surrounding fluid is 
stationary. Three cases, which if fully expanded would 
have jet Mach numbers, Mj , of 1. 1 0, 1.15, and 1 . 1 9, are 
considered. These conditions correspond to the experi- 
mental conditions of Panda [5-8]. In these experiments, 
it was shown that for M } ^ 1.19, the jet noise field is in 
an axisymmetric mode. Hence, the axisymmetric CE/SE 
code described above is appropriate for computation at 
these conditions. 

In the investigation, our attention is focussed on the 
near field of the nozzle since this is the noise source re- 
gion. The diameter, D , of the jet nozzle is chosen as the 
length scale. The density, p 0 , speed of sound, a { j, and 
temperature T 0 in the ambient flow are taken as scales 
for the dependent variables. 

In order to clearly display the upstream propagating 
screech waves, the computational domain was extended 
4 D upstream of the nozzle exit. The full computational 
domain is a circular cylinder of 12 D axial length and 8D 
radius. At the nozzle exit, the inflow plane is recessed by 
two cells so as not to numerically restrict or influence the 
feed-back loop. Ponton et. at .' s [9] straight nozzle lip of 
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NRBC 



Figure 2: Geometry of the computational domain 


0.2D in thickness is adopted in the CE/SE computations. 
Note that in the experimental setup of Panda [5,6], al- 
though the actual nozzle diameter is the same (25.4 mm) 
as in [9], the nozzle exit lip is 0.25 D thick and beveled. 

The unstructured grid currently used is formed by cut- 
ting a rectangle cell into 4 triangles. These rectangu- 
lar cells are non-uniform since good grid resolution is 
needed in the jet shear layer and their numbers in the x 
and y directions are typically 240 and 130, respectively, 
corresponding to about 132000 triangles. The last 10 
cells in the streamwise and radial directions have expo- 
nentially growing x and y size, respectively, and serve as 
buffer zones to essentially eliminate any small remain- 
ing numerical reflection from the outflow boundary. The 
CE/SE scheme is of the a-e type ( 1 3, 1 4 1 with a weighting 
parameter a = 0 and e = 0.5. 

3.1 Initial Conditions 

Initially, the flow of the entire domain is set at the am- 
bient flow conditions, i.e., (using nondimensional vari- 
ables) 


Po — 1* Po — no — 0, vo 0. 

3.2 Boundary Conditions 

At the inlet boundary, the conservative flow variables and 
their spatial derivatives are specified to be those of the 
ambient flow, except at the nozzle exit, where an elevated 
pressure is imposed, i.e ., the jet is under-expanded, as in 
the physical experiments. By using the ideal gas isen- 
tropic relations, it follows that the nondimensional flow 
variables at the nozzle exit, with A/, — 1, are given by 


Pt = 


7(7 + 1 )p c 
2 T r 



2 + (7-1)A//] t 
7 4 - 1 



where T r is the reservoir (plenum) temperature. We will 
also follow the experimental cold-flow condition where 
the reservoir temperature equals the ambient one, i.e., 
Tr = T 

At the symmetry axis, i.e., y — 0, a simple reflective 
boundary condition is applied. At the top and outflow 
boundaries, the Type I and Type II CE/SE non-reflecting 
boundary conditions as described in the next subsection 
are imposed, respectively. The no-slip boundary condi- 
tion is applied on the nozzle walls. 

3.3 Non-Reflecting Boundary Conditions 

In the CE/SE scheme, non-reflecting boundary condi- 
tions (NRBC) are constructed so as to allow fluxes from 
the interior domain of a boundary CE to smoothly exit 
from the domain. There are various implementations 
of the non-reflecting boundary condition and in general 
they have proven to be well suited for aeroaeoustic prob- 
lems [18,19]. The following NRBCs are employed in 
this paper. 

For a grid node (j,n) lying at the outer radius of the 
domain the non-reflective boundary condition (Type I) 
requires that 


(UrYj = (U y ) n j = 0, 

while U r j is kept fixed at the initial steady boundary 
value. At the downstream boundary, where there are 
substantial gradients in the radial direction, the non- 
reflective boundary condition (Type II) requires that 

(u x y; = 0, 


while Uj and (UyYj are now defined by simple extrap- 
olation from the nearest interior node f , i.e., 

u] = u';r l/2 (u y )» = (u y y;r ir2 . 

As will be observed later, these NRBCs, w hen combined 
with the above buffer zones, are robust enough to allow 
a near field computation without disturbing or distorting 
the flow and acoustic fields. 

4 Numerical Results 

In this section, numerical results for the axisymmetric 
underexpanded circular jet are presented and compared 
to experimental results [5-9]. Computations are con- 
ducted for three different jet Mach numbers Mj = LI, 
1.15, and 1.19. With an appropriate time step size, a 
large number (minimum 410,000) of time steps are run 
for these cases in order to achieve appropriate accuracy 
in the Fourier analysis of the time series data. 
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It is important to emphasize that no harmonic forcing 
is imposed in the numerical simulation. The initial im- 
pact of the boundary condition at the nozzle exit stim- 
ulates the jet shear layer and triggers the feedback loop 
that generates the screech waves. 

In the following, the computed jet flow data and acous- 
tic screech data are presented separately, and compar- 
isons are made with available experimental results. 

4.1 Time-Averaged Flow and Shock-Cell 
Structure 

For the typical case of Mj = 1.19, Fig. 3 shows p, p, 
u contours and a numerical Schlieren picture. They are 
time averaged over 40,000 time steps and are obtained 
after the start-up transients have moved out of the com- 
putational domain. In the p, p, and u contours, the shock 
cell structure and the shear layer spreading are clearly 
captured. Note that although the LES model is highly 
simplified, it still distinctly affects the flow field. Since 
the momentum exchange around the shear layer is ac- 
counted for through a high local value of the eddy viscos- 
ity, the shear layer spreads much more rapidly as com- 
pared to our earlier related Euler computations [18], and 
the static pressure maximum increases by 15 — 20%. 

Experimental results for jets are often documented in 
terms of Schlieren pictures. It is straightforward to con- 
struct Schlieren graphs (density-gradient modulus) from 
the numerical results. Such a numerical plot is shown in 
the last panel of Fig. 3. A comparison of the numer- 
ical Schlieren contours from the current computations 
and the experimental Schlieren photograph in [5] (Fig. 4) 
shows good agreement in shock-cell structure. For exam- 
ple, the shock cell width is about 0.8.D in the streamwise 
direction. In both plots, it is observed that the first two 
shock cells appear to be sharp and clear since the shear- 
layer instability wave is too weak at these locations to 
significantly affect the shock cell. However, it interacts 
strongly with the shock cells once it has gained a suffi- 
cient amplitude through its streamwise growth, which is 
evident by the deformation of the third and fourth shock 
cells. The deformation and movement of the shock cells 
are also seen in the experiments [5,6], This is, as pointed 
out by Seiner [2], the mechanism by which the acous- 
tic waves that produce the broadband shock-associated 
noise and the screech tones are generated. Furthermore, 
the fifth shock cell and those further downstream are so 
severely deformed that they almost disappear. 

4.2 Near-Field Radiating Screech Waves 

For the same case of Mj = 1.19, Fig. 5 illustrates the ra- 
diating screech waves at the time level of 410,000 steps. 
Since no forcing is applied in the numerical simulation, 
these waves are a clear indicator of a self-sustained os- 
cillation. For this figure, the very high isobar level con- 
tours, corresponding to hydrodynamic waves around the 


p (.363-. 608) 



Numerical Schlieren 



Figure 3: Time averaged p, p and u contours and numer- 
ical Schlieren plot, Mj = 1.19. 


\i 

k 

Figure 4: Experimental Schlieren plot showing shock 
cell structure. Mj = 1.19. 
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isobars at time step 410,000 



Figure 5: Screech wave radiation, M 3 ~ 1.19. 


jet core area are 'cut off* and the 'colors’ are appropri- 
ately adjusted so that the acoustic waves outside of the 
jet core are clearly displayed. Furthermore, the corre- 
sponding instantaneous numerical Schlieren image is su- 
perimposed on the plot in the core region. The screech 
wavelength is about 1.6D, well in line with experimental 
results (e.g. [5]) and theoretical predictions, e.g. Seiner 
[2], Tam II]. 

4.3 Sound Pressure Level and Screech 
Frequency 

The numerical time history is recorded for a select num- 
ber of locations in the flow field and later post-processed 
to obtain spectral information using Fast Fourier Trans- 
form (FFT) techniques. The data at a selected posi- 
tion is recorded every twentieth actual time step, and the 
recording begins after an initial lime period has elapsed 
(80,000 actual time steps). Typically, each time trace 
contains 33, 000 points and corresponds to more than 0.7 
million computational time steps. Figure 6 is a time trace 
of the unsteady pressure at a location on the nozzle exit 
lip wall (.r = 4.0, y = 0.042) for M 3 = 1.19. The time 
signal is quite irregular. 

Estimates of the sound pressure level (SPL), in dB are 
shown in Figs. 7-9. The abscissas denote frequency in 
Hz. In the following, all the spectral analyses are con- 
ducted at the nozzle exit lip wall (x — 4.0, y — 0.042). 


This is the location used by Ponton et ill [9] in their mea- 
surement, which is somewhat different from the location 
used by Panda [5]. 

Figure 7 displays the SPL for the case of M } = 1.19. 
The SPL shows that there are multiple spikes in the 
spectrum. The two distinct higher spikes correspond to 
screech frequencies of 8890 Hz (SPL= 1 26 dB, .42 mode) 
and 6945 Hz (SPL=131 dB, A\ mode) respectively. 

Figure 8 shows the SPL for the case of A 1 3 = 1.15. 
Still, there are multiple spikes in the spectrum. The 
two distinct higher spikes correspond to screech frequen- 
cies of 9410 Hz (SPL=125 dB, A 2 mode) and 7690 Hz 
(SPL=137 dB, .4] mode) respectively. It is observed 
that these clearly stand out from the mixing and shock- 
associated noise, along with spikes at their first harmon- 
ics. The above two cases support the possibility that 
there could be two screech tone resonances occur at the 
same time or that screech-mode switching is present. 

Figure 9 demonstrates the SPL for the case of AT, = 
1.10. This time, only a single distinct spike stands out 
in the spectrum. The spike corresponds to a screech fre- 
quency of 9040 Hz (SPL=I39 dB, mode). 

4.4 Acoustics Properties Compared to 
Experiments 

Since the feedback loop is very sensitive to any change 
in the nozzle geometry, the computed screech tone fre- 
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Figure 8: SPL at nozzle exit (x = 4.0, y = 0.642), Mj — 
1.15, 75 Hz digital binwidth. 


quency and SPL are slightly different from those ob- 
served in experiments. Even between the experimental 
results of Panda et al. [7] and Ponton et al. [9], there 
are differences in the measured results due to general 
experimental setup and slight differences in nozzle lip 
thickness and shape, although they have the same nozzle 
diameters (25.4 mm) and operating Mach number. 

The numerical results are summarized in Figs. 1 0 and 
1 1, and compared to the experimental data. Figure 10 
displays a comparison of the screech tone frequencies be- 
tween the numerical and experimental results. Both the 
experimental data of Panda et al. [1] and Ponton et al. [9] 
show distinctly the A\ and A 2 screeching modes, and so 
does the numerical data. The CE/SE results agree rea- 
sonably well with Panda et aV s [7] data and seem to 
have a consistent trends for the A\ and A 2 modes. Their 
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difference is less than the difference between the two 
experimental data sets of Panda et al. [7 1 and Ponton 
et al. [9]. Figure 1 1 demonstrates a comparison of the 
screech SPL between the numerical and experimental re- 
sults. The computed SPL results are generally lower than 
the experimental ones, but — except for the Mj = 1.15 
case — are in general agreement with Ponton et al .' s [9] 
data. It is not surprising that the numerical SPL agrees 
better with Ponton et al . ’s [9] data, since their nozzle lip 
(thickness=0.2D, no bevel) is used in the computation 
and the data is taken at the same location ij = 0.642 at 
the nozzle lip wall, while Panda et al.'s [7) data is taken 
at a slightly different location. 

5 Concluding Remarks 

In this paper, the unstructured CE/SE Navier-Stokes 
method is applied for screech tone noise analysis of an 
underexpanded supersonic axisymmetric jet. 

Many hydrodynamic and acoustic aspects of the com- 
puted results are in good agreement with experimental 
findings [5-9], including the variation of screech fre- 
quency with mach number for the A\ and A 2 modes. 
LES modeling was used to better capture the jet spread- 
ing but is not believed to be a critical element for the 
screech feedback loop. 

The advantages of the CE/SE scheme were confirmed 
in the present simulation. The implementation did not 
require any special treatment or parameter selections to 
capture the acoustics. The simple yet robust CE/SE 
NRBC in combination with buffer zones virtually elim- 
inates numerical reflections and allows for an effective 
near-field computation. The use of an unstructured grid 
offers flexibility to deal with more complicated geome- 
tries. 
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